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Influence of strong uniaxial small-scale anisotropy on the stability of inertial- range scaling regimes 
in a model of a passive transverse vector field advected by an incompressible turbulent flow is inves- 
tigated by means of the field theoretic renormalization group. Turbulent fluctuations of the velocity 
field are taken in the form of a Gaussian statistics with zero mean and defined noise with finite corre- 
lations in time. It is shown that stability of the inertial-range scaling regimes in the three-dimensional 
case is not destroyed by anisotropy but the corresponding stability of the two-dimensional system 
can be corrupted by the presence of anisotropy. A borderline dimension d c below which the stability 
of the scaling regime is not present is calculated as a function of anisotropy parameters. 

INTRODUCTION 

It is well known that the breakdown of the classical Kolmogorov-Obuchov phenomenolog- 
ical theory of fully developed turbulence pQ is more noticeable for simpler models of passively 
advected scalar or vector quantity (scalar or vector field) than for the velocity field itself. This 
phenomenon is encoded in the terms intermittency and anomalous scaling . At the same 
time, the problem of a passive advection is easier from theoretical point of view (see, e.g., |3] 
and references therein). Therefore it leads to the fact that the problem of anomalous scaling 
can be understand here in a few ways what is not possible for the present in the problem of 
genuine turbulence. 

One of the most effective approach for studying self-similar scaling behavior is the method 
of the field theoretic renormalization group (RG) 0J|S]. It can be also used in the theory of 
fully developed turbulence and related problems 0017], e.g., in the problem of a passive 
scalar (or vector) field advected by a given stochastic environment. 

In |Hj the field theoretic RG was applied to the so-called rapid-change model of a passive 
scalar advected by a self-similar white-in-time velocity field which is also known as Kraichnan 
model. It was shown that within the field theoretic RG approach the anomalous scaling is 
related to the existence of " dangerous" composite operators with negative critical dimensions 
in the framework of the operator product expansion (OPE) [SJ|H1[7]. 

Afterwards, various generalized descendants of the Kraichnan model, namely, models with 
inclusion of large and small scale anisotropy, compressibility, and finite correlation time of 
the velocity field were studied by the field theoretic approach (see [5] and references therein). 
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Moreover, advection of a passive vector field by the Gaussian self-similar velocity field (with 
and without large and small scale anisotropy, pressure, compressibility, and finite correlation 
time) has been also investigated and all possible asymptotic scaling regimes and crossovers 
among them have been classified |l(JllllllT2lll3) . General conclusion is: the anomalous scaling, 
which is the most important feature of the Kraichnan rapid-change model, remains valid for 
all generalized models. 

In what follows we shall begin with investigation of one particular model of a passive vector 
advected by a Gaussian velocity field with finite correlation time in the presence of the small- 
scale anisotropy, namely the model where the stretching term is absent (the so-called A = 
model, see, e.g, [111 113) ~). This model is specific from several points of view but maybe the 
most important fact is that in contrast to the other models of passive vector admixture where 
the anomalous scaling is related to the composite operators built of the vector field without 
derivatives |121 ITS] in the case under consideration it is related to the composite operators 
built solely of the gradients of the field. This fact radically changes the complexity of the 
problem especially in the anisotropic case (see, e.g., |11M14| and references therein). Thus, in 
some sense, it can be consider as a further step to the nonlinear Navier-Stokes equation. In 
what follows, we shall present only the first part of the RG analysis, namely, we shall analyze 
the influence of the small-scale anisotropy on the infrared (IR) stability of the possible scaling 
regimes of the model. It will be seen that complexity of this task is also very close to the 
corresponding problem in the stochastic Navier-Stokes equation 15 . 

1. FIELD THEORETIC FORMULATION OF THE MODEL 

We shall consider the model of the advection of transverse (solenoidal) passive vector field 
b = b(x, t) which is described by the following stochastic equation 

d t b = u Ab - (v • V)b + f , (1) 

where dt = d/dt, A = V 2 is the Laplace operator, vq is the diffusivity (a subscript denotes 
bare parameters of unrenormalized theory), and v = v(x, i) is incompressible advecting ve- 
locity field. The vector field f = f (x, t) is a transverse Gaussian random (stirring) force with 
zero mean and covariance 

D{ j = (f i (xL,t)f j (x',t'))=5(t-t')C ij (r/L), r = x - x' (2) 

where parentheses (...) hereafter denote average over corresponding statistical ensemble. The 
noise defined in Eq. J5J maintains the steady-state of the system but the concrete form of the 
correlator will not be essential in what follows. The only condition which must be fulfilled by 
the function CV,(r/L) is that it must decrease rapidly for r = |r| L, where L denotes an 
integral scale related to the stirring. 

In real problems the velocity field v(x) satisfies Navier-Stokes equation but, in what follows, 
we shall work with a simplified model where we suppose that statistics of the velocity field is 
given in the form of Gaussian distribution with zero mean and pair correlation function 

f H^]<rlt ) 

(viix^x'V^D^x') = J ^(kFKkK*')^), (3) 



3 



where d is the dimension of the space, k is the wave vector, and Rij (k) is a transverse projector. 
In our uniaxial anisotropic case it is taken as (see, e.g., |12j and references therein) 

Rij (k) = (1 + a x (n • k) 2 /fc 2 ) Pij (k) + a 2 n sni P ts (k)P j7 (k) , (4) 

where Py(k) = 5ij — kikj/k 2 is common isotropic transverse projector, the unit vector n de- 
termines the distinguished direction, and ai, ai are parameters characterizing the anisotropy. 
From the positiveness of the correlation tensor one immediately finds restrictions on the 
values of the above parameters, namely ai,2 > — 1. The function D v (uj,\s) in is taken in 
the following form [T3"| 

n ii r/3 1.4— d— 2s — n 

D^k) = - ^ 90U h Zt T 5^, (5) 

where go plays the role of the coupling constant of the model (a formal small parameter 
of the ordinary perturbation theory), the parameter uo gives the ratio of turnover time of 
scalar field and velocity correlation time, and the positive exponents e and r\ are small RG 
expansion parameters. The coupling constant go and the exponent e control the behavior of 
the equal-time pair correlation function of velocity field and the parameter uq together with 
the second exponent r] are related to the frequency to ~ ito^ofc 2- '' which characterizes the 
mode k. The value e = 4/3 corresponds to the celebrated Kolmogorov "two-thirds law" for 
the spatial statistics of velocity field, and 77 = 4/3 corresponds to the Kolmogorov frequency. 
Simple dimensional analysis shows that go and uq, which we commonly term as charges, are 
related to the characteristic ultraviolet (UV) momentum scale A (or inner legth I ~ A -1 ) by 

. 90 ~A 2£ , uo^A". (6) 

The stochastic problem can be treated as a field theory with action functional 0J|S] 

5($) = bj [(-d t - Vidi + vo A + i^oXio(n ■ d) 2 ) 5 jk + n 3 u (X20A + x 3 o(n ■ d) 2 ) n k ] b k 

— (^[^■rs-^aj), (?) 

where and d(- are given in © and respectively, b' is an auxiliary vector field (see, 
e.g., 0), and the required integrations over x = (x, t) and summations over the vector indices 
are implied. In action J7J) the terms with new parameters XiO;X20> and X30 are related to the 
presence of small-scale anisotropy and they are necessary to make the model multiplicatively 
renormalizable. Model (JZJ) corresponds to a standard Feynman diagrammatic technique (see, 
e -g-j jl2M15j for details) and the standard analysis of canonical dimensions then shows which 
one-irreducible Green functions can possess UV superficial divergences. 

The functional formulation J7| gives possibility to use the field-theoretic methods, includ- 
ing the RG technique to solve the problem. By means of the RG approach it is possible 
to extract large-scale asymptotic behavior of the correlation functions after an appropriate 
renormalization procedure which is needed to remove UV-divergences. 
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2. SCALING REGIMES OF THE MODEL 



5 

T|4 

3 
2 
1 

-1 
-2 
-3 
-4 
-5 



FPI 


2e=r,/ FPU / 


FPIII 


FPIV 



-5 -4 -3 -2 -1 1 2 3 4 5 



Fig. 1: The scaling regimes of the model in the 
e — r\ plane. The region FPI corresponds to the 
trivial rapid-change limit: g*/u* = 0, 1/u* — 0. 
The region FPU corresponds to the the non-trivial 
rapid change limit: g*/u* > 0,1/it* = 0. The re- 
gion FPIII corresponds the trivial "frozen" limit: 
g* = 0, u* = 0. The region FPIV is related to the 
non-trivial "frozen" limit: g* > 0, it* = 0. In the 
end, the line e — rj (FPV) corresponds to the more 
interesting scaling regime with g* > 0, < it* < oo. 



Details of the one-loop field-theoretic RG 
analysis of the model will be done elsewhere. 
Here we only conclude that using the RG 
analysis leads to the following result: possi- 
ble scaling regimes are given by the IR sta- 
ble fixed points of the system of five non- 
linear RG differential equations (flow equa- 
tions or also known as Gell-Mann-Low equa- 
tions) for five scale dependent effective vari- 
ables (charges) C = {g, u, xi, x.2, X3j of the 
model which are functions of the dimension- 
less scale parameter t = k/ A [5J ^] . The 
system of the flow equations is defined by 
the so-called /3-functions of the model (they 
are functions of the charges, anisotropy pa- 
rameters, space dimension, and parameters 
e, rj) and it has the following form 



dt 
du 



t 



dt 
dt 



g = g(-2e + 2 7l ) , (8) 

/3 u =n(-77 + 7i), (9) 

^*=Xi(7i-7i+i) ! (10) 
i= 1,2,3, 



where the functions ji, i = 1, 2, 3, 4 are given 
by the following expressions (the one-loop 
approximation) 



7i 
7*+i 
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(11) 
(12) 



where Sd = 2n d / 2 /T(d/2) is the surface of the d dimensional sphere, w\ = (1 + u + xix 2 ), u>2 = 
(l + u + xix 2 + (X2 +X3^ 2 )(1 - x 2 )) 7 and the coefficients Ki(x 2 , g, u, Xi, X2, X3, «i, "2, d), i 
1 . 2, 3, 4 are some huge polynomials in respect to all variables and their explicit form will be 
given elsewhere. The scale parameter t belongs to the interval < t < 1 with the initial 
conditions given at t = 1 and the IR stable fixed point corresponds to the limit t — > 0, i.e., 
C\ t=0 = C*. 

We have performed a numerical analysis of this system of differential equations and our aim 
was twofold. First of all, we have found all possible scaling regimes and we have analyzed the 
regions of their IR stability in the e — rj plane. The results of this analysis are shown in Fig.^ 
where it is shown that the model exhibits five different scaling regimes (two for rapid-change 
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Fig. 2: Dependence of the borderline dimension d c on the parameters ai and ai for different values 
of u = u* . The corresponding scaling regime is stable above the given curve. 

limit, two for so-called "frozen" limit, and one general with nonzero u») (see, e.g., and 
references therein) . The second question which was investigated is related to the dependence 
of stability of the above mentioned scaling regimes on the anisotropy parameters ol\, ol^ and 
on the dimension of the space d. We have found the so-called borderline dimension d c between 
stable and unstable regimes as a function of anisotropy parameters ai,a2 and parameter u*. 
The results are shown in Fig.[21 for some special situations. One can see that the presence 
of small-scale anisotropy leads to the violation of the stability of the corresponding scaling 
regimes below d c £ [2, 3] for appropriate values of anisotropy parameters. But from the point 
of view of further investigation of anomalous scaling the most important conclusion is that 
all the three-dimensional scaling regimes remain stable under influence of small-scale uniaxial 
anisotropy. 

CONCLUSIONS 

Using the field theoretic RG we have studied the influence of small-scale uniaxial anisotropy 
on the stability of the scaling regimes in the model of a passive vector advected by given 
stochastic environment with finite time correlations. It is shown that the system exhibits five 
different scaling regimes. They are related to the values of the parameters e and rj. On the 
other hand, the stability of all these scaling regimes are influenced by presence of small-scale 
anisotropy which is demonstrated in the existence of the so-called borderline dimension d c 
which is a function of the anisotropy parameters. The d c is defined as dimension above which 
the corresponding scaling regime is stable and below which the stability of the regime is de- 
stroyed. All calculations have been done at the one-loop level. The results will be used in the 
further investigations of the anomalous scaling of the model. 
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